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ABSTRACT 

This  paper  presents  a  computer  program  for  a  common  type  of  analysis 
of  the  slope  stability  problem,  viz. ,  the  possibility  of  slope  failure  by 
translation  of  a  massive  block  along  a  weak  layer  of  soil.  The  problem, 
which  can  occur  in  either  natural  or  man-made  slopes ,  is  most  generally 
referred  to  as  the  "sliding  block  problem". 

Variation  in  the  water  surface  position  requires  three  subroutines 
or  cases.  The  program  automatically  sequences  selected  potential  sliding 
surfaces  one  by  one,  then  selects  the  desired  water  surface  case,  and 
finally  computes  the  factor  of  safety  against  sliding  along  the  base  of 
the  central  block. 

The  analysis  is  based  on  total  unit  weights  and  boundary  forces. 
It  is  possible  to  consider  ten  different  soil  types  having  very  different 
soil  parameters,  viz.,  unit  weight,  Mohr-Coulomb  cohesion  intercept  and 
Mohr-Coulomb  angle  of  friction.  A  maximum  of  twelve  continuous  soil  lay- 
ers at  any  inclination  can  be  considered  in  the  present  program.  A  total 
of  ten  vertical  strip  loads  of  different  intensities  can  be  placed  on  the 
ground  surface  anywhere  below  the  toe  and  above  the  crest.  Finally,  with 
all  the  above  information,  ten  sliding  surfaces  can  be  concurrently  ana- 
lyzed for  the  factor  of  safety.  Said  factor  is  applied  to  the  strength 
of  the  soil  at  the  base  of  the  central  block,  assuming  limiting  equilibrium 
for  the  active  and  passive  earth  pressure  forces  at  the  ends  of  the  central 
block . 

The  paper  is  complete  with  a  flow  chart  of  the  computer  program  and 
two  illustration  problems. 


)han,  A. 


LIST  OF  FIGURES 

Figure  Page 

1.  Slope  in  Stratified  Soil  Profile  2 

2.  General  Problem  of  Sliding  Block  with 

Submerged  Water  (Case  2)  h 

3.  Analysis  of  Forces  on  Active  Wedge  (Case  2)   7 

U.  Analysis  of  Forces  on  Passive  Wedge  (Case  2) 9 

5.  Forces  on  Sliding  Block  (Case  2) 11 

6.  Flow  Chart  of  Block  Program lU 

7.  Illustration  Problem  No.  1   16 

8.  Illustration  Problem  No.  2   18 


INTRODUCTION 

The  stability  of  man-made  and  natural  slopes  has  always  been  an  im- 
portant topic  of  discussion  in  the  field  of  civil  engineering.   Yet,  fail- 
ure of  man-made  fills  and  cuts  probably  occurs  more  frequently  than  all 
other  failures  of  civil  engineering  structures  combined.   Although  an 
understanding  of  the  major  factors  which  contribute  to  failure  of  slopes 
has  improved  considerably,  our  predictive  ability  remains  less  than 
satisfactory. 

This  paper  addresses  the  problem  of  the  "sliding  block",  i.e.,  an 
essentially  rigid  mass  sliding  in  a  weak  layer.  At  first  glance,  this 
seems  to  be  a  rather  simple  problem.   But  when  practical  variations  in 
soil  profile  are  considered,  as  well  as  water  levels,  boundary  geometries 
and  loadings,  and  uncertainties  of  position  and  shape  of  the  most  critical 
sliding  surface,  the  solutions  require  reasonably  large  computer  systems. 

When  a  slope  is  underlain  by  one  or  more  strata  of  very  soft  or 
loose  materials ,  the  most  critical  sliding  surface  may  not  be  even  approx- 
imately circular,  as  shown  in  Figure  1.  Rather  there  is  a  3-plane  surface 
of  potential  sliding  in  which  a  maximum  amount  of  the  surface  lies  within 
the  weak  material. 

An  initial  programmed   solution  by  one  of  the  junior  authors 
(Mendez  (1972))  was  quite  general  with  respect  to  the  shape  of  the  3-plane 
surface,  but  to  accommodate  this  feature  the  profile  was  simplified  to 
two  soil  layers,  viz.,  a  strong  soil  over  a  weak  one.  A  second  program, 
reported  in  this  paper,  makes  simplifying  assumptions  with  respect  to 
the  shape  of  the  sliding  surface,  but  is  quite  versatile  with  respect  to 
the  profile  and  boundaries.   This  second  program  seems  to  better  meet  the 
analytical  requirements  of  the  Indiana  State  Highway  Commission. 
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Location  of  failure  surface  depends  on 
relative  strength  and  orientation 
of  layers. 


Stable  soil; 
high  strength 


FIGURE  I 


SLOPE  FAILURE  BY  SLIDING 

The  type  of  failure  usually  assumed  in  slope  stability  analysis  is 
the  one  piece  slide  (HRB  (1958))  .  The  failure  is  one  in  which  the  moving 
body  is  essentially  rigid  and  the  failing  mass  is  separated  from  the 
unmoved  one  by  a  surface  of  assumed  shape.  Where  the  soil  is  grossly 
homogeneous  it  seems  logical  that  the  failure  surface  would  be  roughly 
circular,  and  in  the  interest  of  simplicity  it  is  usually  made  exactly 
so.  A  recent  overview  of  the  circular  type  analysis,  by  the  well  known 
methods  of  slices,  is  contained  in  Carter,  Lovell  and  Harr  (1971). 

Where  there  is  evidence  of  definite  differences  in  shearing  resistance 
in  the  soil  profile,  it  is  well  to  consider  potential  failure  surfaces 
which  follow  the  surfaces  of  weakness.  Several  methods  of  handling  irreg- 
ular surfaces  are  reported  by  Morganstern  and  Price  (1965),  Carter,  Lovell 
and  Harr  (1971),  and  Mohan  (1971). 

A  special  case  of  the  irregular  sliding  surface  has  been  shown  in 
Figure  1,  where  the  potential  failure  planes  have  a  maximum  length  in  the 
weaker  materials.  The  potential  failing  block  is  actually  a  combination 
of  active  and  passive  wedges,  with  a  central  trapezoidal  block  based  in  a 
weak  layer.  Examples  of  simplified  solutions  to  this  problem  are  given 
in  Department  of  the  Navy  (1971)  and  United  States  Steel  (1972),  as  well  as 
Mendez  (1972). 

A  GENERAL  SOLUTION  TO  THE  SLIDING 
BLOCK  PROBLEM 

Figure  2  shows  the  free  body  diagram  with  a  full  quota  of  complexities 
in  boundary  geometries  and  forces,  i.e.,  these  could  be  more  simple  in  a 
given  instance.  Incorporation  of  a  water  surface  and  associated  water 
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forces  into  the  problem  makes  it  convenient  to  consider  three  cases ,  each 
with  its  appropriate  subroutine  in  the  computer  solution.  The  upper 
boundary  slopes  reading  left  to  right  in  Figure  2  are  referred  to  as  the 
"down  slope",  the  "middle  slope"  or  simply  the  "slope"  and  the  "upper 
slope".  The  cases  are: 

Case  1.  When  the  water  surface  is  below  the  trial  sliding  surface. 

Case  2.  When  the  water  surface  is  partly  above  and  partly  below  the 
ground  surface,  but  above  the  trial  sliding  surface. 

Case  3.  When  the  water  surface  is  anywhere  below  the  ground  surface, 
but  above  the  trial  sliding  surface. 

It  is  assumed  that  the  right  hand  wedges  are  in  a  state  of  limiting 
active  earth  pressure  and  the  left  hand  wedges  are  in  a  state  of  limiting 
passive  earth  pressures.  Simplifying  assumptions  are  employed  with  re- 
spect to  the  inclinations  of  the  wedges  surfaces  and  the  directions  of 
the  earth  pressure  forces.  Although  the  right-hand  and  left-hand  wedges 
are  assumed  to  be  on  the  verge  of  sliding,  there  is  in  general,  an  incom- 
plete mobilization  of  the  shearing  resistance  along  the  base  of  the  block, 
i.e.,  the  factor  of  safety  is  defined  with  respect  to  the  shearing 
resistance-shearing  force  ratio  along  this  surface. 

The  wedge  inclination  and  earth  pressure  force  direction  assumptions 
are  those  which  apply  for  a  simple  Rankine  case.  They  are  employed  by 
others  (Dept.  of  the  Navy  (l9fl)),  and  have  been  shown  to  be  good  approxi- 
mations of  the  most  critical  values,  for  a  number  of  cases  tested  by 
Mendes  (1972). 

To  be  certain  that  all  assumptions     inherent  to  the  solution  are 
understood,  they  are  listed  in  detail  below. 


1.  Problem  is  two  dimensional. 

2.  The  ground  surface  is  defined  by  three  slopes,  and  a  well  defined  toe 
and  crest. 

3.  Soil  strata  are  laterally  continuous. 

h.     Soil  properties  in  layers  are  defined  by  y,  c,  and  <|>  (where  c  or  I 
can  be  equal  to  zero). 

5.  Sliding  surface  at  the  base  of  the  block  and  between  the  slide  wedges 
is  a  plane. 

6.  All  lateral  forces  on  vertical  wedge  boundaries  are  normal  to  these 
boundaries,  i.e.,  there  are  no  shear  forces  on  these  boundaries. 

7.  The  factor  of  safety  is  figured  for  the  base  of  the  sliding  block  only. 
The  movement  required  to  mobilize  limiting  active  and  passive  pressures 
is  smaller  than  the  movement  required  to  mobilize  the  shearing  strength 
of  the  weak  soil  strata. 

8.  The  wedge  slip  surfaces  are  at  (1*5  +  <j>/2)  and  (1*5  -  4>/2)  with  the  hori- 
zontal for  active  and  passive  wedges,  respectively. 

9.  The  active  and  passive  forces  are  computed  by  satisfying  static  equi- 
librium and  after  assumptions  6  and  8. 

10.   Seepage,  if  any,  is  in  a  steady  state.   However,  water  pressures  are 
calculated  at  any  point  as  if  they  were  hydrostatic. 

The  analysis  of  forces  is  demonstrated  in  Figures  2,  3,  and  h.     The 
analysis  is  divided  into  three  parts : 

(i)  calculation  of  forces  on  central  block  due  to  active  wedge; 
(ii)  calculation  of  forces  on  central  block  due  to  passive  wedge;  and 
(iii)  calculation  of  base  forces  on  central  block  and  of  the  factor  of 
safety  against  sliding  along  this  base. 
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The  analysis  of  forces  is  illustrated  for  water  surface  case  2,  but 
the  other  cases  follow  directly. 

Figure  2  shows  a  rather  complex  problem  space  section,  with  multiple 
soil  layers  at  variable  inclinations  and  with  very  different  soil 
properties . 

(i)  Analysis  of  Active  Forces  on  Central  Block 

Figure  3  shows  the  active  wedge  from  Figure  2,  divided  into  small 
wedges  governed  by  the  intersection  of  the  assumed  slip  surface  and  soil 
boundaries . 

Consider  a  typical  polygon  of  forces  for  any  (nth)  wedge  in  Figure  3. 

Summation  of  all  the  forces  in  the  x  and  y  directions  and  equating  to 

zero  yields  the  following  equations , 

ZF  =  0 
x 

PAn  ■  UARn  -  UALn  -  UAnCos   (1+5  -  <J>n/2)  +  NA'n  Cos    (1+5  +  <f>n/2) 

-  CAn  Cos   (U5  +  <l»n/2)  (l) 

ZF     =  0 

y 

WAn  «  CAn  Sin  (1+5  +  <J>n/2)  +  UAn  Sin   (1+5  *  <frn/2) 

+  NA'n  Stn  (1*5  +  <j>n/2)  (2) 

Elimination  of  NA'n  from  equations   (l)  and  (2)  yields  an  expression 
for  the  incremental  active  force  for  the  nth  wedge, 
PAn  =  WAn  Tan   (1+5  -  <j>n/2)  -  2  CAn  Cos   (1+5  +  <jm/2) 

+   (UARn  -  UALn)  +  UAn    |cos   (1*5  -  <(>n/2) 

-  Tan  (U5  -  <frn/2)  Sin  (1+5  -  <frn/2)}  (3) 
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(ii)  Analysis  of  Passive  Forces  on  Central  Block 
Figure  k   shows  the  forces  acting  on  the  passive  wedge  from  Figure 
2.   Consider  a  typical  polygon  of  forces  acting  for  an  nth  passive  wedge 
in  Figure  1*.  Sum  forces  in  the  x  and  y  directions,  and  equate  to  zero 
for  equilibrium. 

EF  =  0 
x 

PPn  ■  UB     Sin  6,  +  CP     Cos   (1*5  -  <t>n/2)  +  ULn  +  UPn  Cos    (1*5  +  <*V2) 
n  J.  n 

-  URn  +  NP'n    Cos  (1*5  -  <jm/2)  CO 

EF     ■  0 

y 

WPn  =  NP'n  Sin  (1*5  -  <|>n/2)  -  U3Q  Cos  $x  -  CPn  Sin  (1*5  -  4>n/2) 

+  UPn  Sin   (1*5  +  <fV2)  (5) 

Elimination  of  NP'n  from  equations   (h)  and  (5)  yields  an  expression 
for  the  incremental  passive  force  for  the  nth  wedge. 

PPn     «  WPn  Tan   (1*5  +  <Jm/2)  +  2  CPn  Cos   (1*5  -  4>n/2) 

+  Ug     |  Sin  31  +  Cos  &±  Tan   (1*5  +  <j>n/2)^ 

+   (ULn  -  URn)  +  UPn  jcos    (1*5  +  <J>n/2)  -  Sin   (1*5  +  <J>n/2)  Tan   (1+5  +  <jm/2)| 

(6) 

(iii)  Analysis  of  Forces  on  Central  Block  and  Calculation  of  Factor 
of  Safety 

Figure  5  shows  the  appropriate  free  body  from  Figure  2.  The  factor 
(FS)  is  commonly  called  the  factor  of  safety,  although  it  is  better  in- 
terpreted as  a  strength  reduction  factor,  i.e.,  if  the  real  strength  were 
divided  by  this  factor,  a  reduced  strength  would  obtain  at  which  failure 
would  impend.  Note  that  the  base  sliding  surface  can  be  inclined  up 
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(8   )   or  down  (9   )  with  respect  to  the  horizontal,  or  may  be  horizontal 
(9  =  0). 

For  9     and  where  forces  are  summed  normal   (N)    and  tangential   (9)   to 
the  sliding  surface, 

ZFN=° 

NB'  +  UBP  ■  PAA  Sin  9  -  PPP  Sin  9  +  WB  Cos  9  +  UBH  Cos  &2  Cos  9 

-  UBL  Sin  9  +  UBR  Sin  9  -  UBH  Sin  62  Sin  9  (7) 

ZFe  =  0 

g  +  NB'  Tan  0  =  pAA  Cos  Q  _  ppp  Cos  Q  _  v^  Sln  0 

-  UBH  Cos  $2  Sin  9  -  UBH  Sin  8  Cos  9 

-  UBL  Cos  9  +  UBR  Cos  9  (8) 

Elimination  of  NB'  from  equations  (7)  and  (8)  yields  an  expression 
for  the  factor  of  safety  for  a  particular  trial  sliding  surface, 


CB  +  (PAA  Sin  9  -"PPP  Sin  9  +  WB  Cos  9  +  UBH  Cos  B2  Cos  9 


(9) 


-  UBL  Sin  9  +  UBR  Sin  9  -  UBP  -  UBH  Sin  $2  Sin  8)  Tan  <f> 
^  =  (PAA-PPP)  Cos  9  -  WB  Sin  9  -  UBH  Cos  62  Sin  9 

-  UBH  Sin  32  Cos  9  -  UBL  Cos  9  +  UBR  Cos  9 

For  9+ 

CB  +  (PPP  Sin  9  -  PAA  Sin  9  +  WB  Cos  9  +  UBH  Cos  82  Cos  9 

+  UBL  Sin  8  -  UBR  Sin  9  -  UBP  +  UBH  Sin  62  Sin  9)  Tan  <f> 
re  ~ (PAA  -  PPP)  Cos  9  +  WB  Sin  9  +  UBH  Cos  62  Sin  9  (l0) 

-  UBH  Sin  32  Cos  9  -  UBL  Cos  9  +  UBR  Cos  9 

For  9=0  (Horizontal  slope) 

CB  ■»■  (WB  -  UBP  +  UBH  Cos   Bg)   Tan  <fr 
re  '   (PAA-PPP)   -  UBH  Sin  0„  -  UBL  +  UBR  (l1* 
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THE  COMPUTER  PROGRAM  AND  ITS  CAPABILITIES 

The  flow  chart  for  the  program  is  shown  in  Figure  6.   The  program 
has  been  written  in  FORTRAN  IV  language,  and  at  present  it  is  workable  on 
the  CDC  65OO  computer.   It  is  made  up  of  a  main  program  and  six  supporting  sub- 
routines. The  program  makes  use  of  common  storage  to  optimize  usage  of 
high  speed  core  and  minimum  computation  time. 

The  program  is  capable  of  handling  the  following  variables, 
(i)  Multiple  (up  to  11)  continuous  soil  layers  at  any  inclination; 
layer  boundaries  are  straight. 

(ii)  Top  ground  surface  made  up  of  three  slopes  and  well  defined  toe 
and  crest  points. 

(iii)  Soil  properties  defined  by  y,  c  and  4>  (c  or  <f>  can  equal  zero), 
(iv)  Multiple  (up  to  10)  uniform  strip  loads  on  ground  surface  of  the 
upper  and/or  lower  slopes. 

(v)  Water  surface  anywhere  in  the  problem  space.   The  water  surface 
is  defined  by  continuous  straight  lines  and/or  by  a  non-linear  surface 
defined  by  7  or  less  known  coordinates. 

(vi)  Multiple  trial  sliding  surfaces  at  the  bottom  of  the  central 
block.  These  can  be  at  any  inclination  and  as  many  as  10  can  be  analyzed 
in  a  single  run. 

Specific  trial  surfaces  are  input  for  analysis.  No  searching  tech- 
nique (for  identification  of  a  minimum  FS)  is  recommended,  although  some 
ideas  on  this  are  contained  in  Carter,  Lovell  and  Harr  (1971). 

The  active  and  passive  force  subroutines  are  potentially  valuable  of 
themselves  in  the  solution  of  lateral  earth  pressure  problems. 
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ILLUSTRATION  PROBLEMS 

The  purpose  of  the  illustration  problems  is  three  fold: 
(l)  to  demonstrate  the  use  of  the  computer  program;  (2)  to  show  the 
versatility  and  several  options  of  the  program;  and  (3)  to  serve  as 
a  check  for  duplicated  decks.  Two  separate  hypothetical  problems  are 
chosen  for  this  purpose. 

Illustration  Problem  No.  1 

This  first  illustration  problem  involves  a  simple  soil  profile 
shown  in  Figure  7.  Solutions  are  obtained  for  three  central  block 
sliding  surfaces  and  for  three  location  of  the  water  surface  for 
each  sliding  surface.  The  results  are  given  in  Table  1. 

Illustration  Problem  No.  2 

This  second  problem  is  more  complex  and  is  shown  in  Figure  8.  This 
problem  is  also  solved  for  three  slopes  of  sliding  surfaces  in  combination 
with  three  locations  of  the  water  surface  for  each  sliding  surface. 
The  results  are  again  shown  in  Table  1. 
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TABLE  1 


Summary  of  Results  for  Two  Illustration  Problems 


Case 
Analysed 


Number  and 
Slope  of 
Sliding  Surface 


Factor  of  Safety 


Problem  No.   1 


Problem  No.    2 


Case  1 


Case  2 


Case  3 


1. 

et 

2. 

0  =  0 

3. 

e" 

1. 

e+ 

2. 

8  =  0 

3. 

e" 

1. 

e+ 

2. 

6  =  0 

3. 

e" 

1.91* 
1.87 
1.83 

1.66 
1.65 
1.6U 

1.81* 
1.52 
1.66 


2.07 
2.2U 
2.U2 

1.83 
1.98 
2.12 

1.81+ 

1.8l 
1.75 
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CONCLUSIONS  AND  RECOMMENDATIONS 

The  primary  objective  of  this  research  was  the  development  of  a 
computer-assisted  system  for  rapid  prediction  of  the  factor  of  safety 
of  slopes  where  the  mode  of  failure  is  a  sliding  block.  The  resulting 
program  is  sufficiently  versatile  to  accommodate  a  3-slope  ground 
surface  and  a  subsurface  profile  with  spatial  variations  in  material 
properties,  a  steady  state  flow  domain,  and  uniform  strip  ground 
surface  loadings.  Up  to  ten  trial  sliding  surfaces  can  be  analyzed 
concurrently,  with  the  base  of  the  central  wedge  at  any  inclination 
in  any  selected  soil  layer.  The  program  automatically  sequences 
the  trial  sliding  surfaces,  computing  a  factor  of  safety  for  each. 
Since  many  sliding  surfaces  may  have  to  examined,  i.e.,  there  is  no 
systematic  search  technique  which  assures  identification  of  a  minimum, 
this  is  a  most  important  feature. 

It  was  desired  to  develop  a  system  which  chould  be  used  on 
smaller  computers.   Consequently,  the  program  uses  a  small  storage 
and  short  computation  time. 

Hopefully,  the  program  will  enable  a  designer  to  check  against 
this  mode  of  instability  for  any  slope  where  there  is  reason  to 
suspect  that  it  may  occur.  Such  suspicision  would  ordinarly  accrue 
from  study  of  boring  logs  and  profiles.  Sliding  blocks  can  be 
based  in  any  soil  stratum  of  below-average  strength.  When  there  is  no 
evidence  of  weak  layers,  it  is  likely  that  some  common  form  of  the 
circular  or  rotational  slump  analysis  will  be  employed.  In  questionable 
cases,  both  types  of  analysis  may  be  undertaken  and  factors  of  safety 
compared. 
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Any  computer  program  should  be  tested  for  reliability  by  generation 
of  solutions  to  common  problems  through  different  programs  or  manual 
calculations.   Unfortunately,  this  is  usually  possible  for  only  simple 
examples,  since  the  motivation  for  development  of  the  new  programs  is 
an  inability  of  old  ones  to  accommodate  the  desired  level  of  physical 
complexity. 
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